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algorithms  are  presented  for  both  Burmester  and  Conic 
Theory. 

An  extension  of  Burmester  Theory  to  the  synthesis 
of  six  point-positions  of  the  moving  plane  is  presented 
with  a theorem  of  uniqueness  for  six  positions  of  the 
moving  plane.  The  nine  precision  point  problem  is 
reduced  to  a solution  of  five  equations  with  five  unknown 
variables. 

Tabulated  coefficients  for  the  analysis  and  synthesis 
algorithms  provide  a systematic  and  simplified  means  of 
computation  for  all  case  studies. 


CHAPTER 


GENERAL  BACKGROUND 

In  1967  the  author  completed  research  which  coalesced 
the  study  of  finite  (positional)  and  infinitesimal 
(instantaneous)  kinematics  by  presenting  an  algorithm  for 
five  multiply  separated  positions  (any  multiple  combination 
of  finite  and/or  infinitesimal  positions)  in  coplanar 
motion.  The  work  is  similar  to  that  outlined  by  Bottcma 
(1)  for  the  finite  case  and  requires  the  resolution  of 
four  linear  and  two  non-linear  algebraic  equations  based 
on  circular  constraints  in  the  fixed  plane.  A generalised 
computer  program  capable  of  solving  all  seven  cases  was 
the  culmination  of  this  work.  Since  its  completion  there 
have  been  few  research  publications  associated  with 
generalised  positional  synthesis  for  Burmester  or  Conic 
constraints.  Consequently,  as  in  many  disciplines  of 
study,  kinematics  appeared  to  have  reached  a plateau  in 
its  capabilities  in  coplanar  synthesis.  Often  when  this 
happens  the  only  format  to  follow  is  to  survey  the  field 
of  discipline  with  new  conceptual  ideas  to  see  if  this 
might  precipitate  different  concepts  which  will  reflect 
upon  ar.d  advance  the  particular  area  of  study.  Success  is 


primarily  dependent  upon  taking  an  interrogative  viewpoint 
in  the  light  of  what  is  known  and  what  has  created  such  a 
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TYPES  OF  CONICS 

1)  ELLIPSE 

2)  PARABOLA 

3)  HYPERBOLA 

4)  CIRCLE  (SPECIAL  ELLIPSE) 


Figure  (1-J 


Analysis  of  the  Coupler  Curve  Function 
for  Conic  Constraint  Pairs’ 

In  recent  years  many  kinematic  researchers  have 
expended  considerable  effort  to  define  and  analyze  the 
output  function  (coupler  curve)  of  the  simple  four-bar 
mechanism.  Since  this  study  is  concerned  with  conic 
constraint  pairs  in  motion,  an  algorithm  is  presented  to 
coalesce  and  generalize  the  analysis  of  the  conic 
mechanisms  with  the  mode  of  synthesis. 

As  will  be  shown  in  the  text,  the  analysis  of  the 
circular  constraint  system  employs  the  same  algorithm  as 
the  generalized  conic  constraint  system.  However,  it 
becomes  advantageous  to  discuss  the  two  systems  separately 
because  of  the  difference  in  the  complexity  of  the  solutions. 

The  basic  concern  for  the  develop  rent  of  an  analysis 
algorithm  was  to  acquire  a method  of  verifying  the  results 
provided  by  the  synthesis  program.  The  author  found  that 
his  efforts  were  markedly  reduced  as  a result  of  this 
particular  investigation  and  formulation. 

Analysis  of  the  Coupler  Curve  Function 
for  Burmester  Constraints 

from  a kinematic  standpoint,  analysis  of  coplanar 
linkages  is  cerfainly  a well-known  area  of  study.  Works 
by  Freudenstein  [21  have  provided  explicit  results  for  all 
types  uf  circular  (Burmester)  constrained  systems.  There- 
fore, it  is  not  the  intention  of  this  paper  to  initiate  new 
areas  of  research  for  analysis  but  simply  to  develop  a 


discipline  or  format  which  is  compatible  with  the  described 
mode  of  synthesis.  Methods  will  be  discussed  for  analysing 
positions  as  well  as  orders  of  contact  (tangent,  curvature, 
etc.)  for  Burme3ter  constraints. 

Linkage  Synthesis  for 
Burmester  and  Conic  Constraint  Pairs 

Pecent  works  by  Sparks  and  Tesar  (3)  have  provided  a 
mode  of  synthesizing  a Burmester  pair  while  works  by 
Hain  [ 4 J have  provided  a graphical  mode  of  synthesizing 
the  ground  link  for  Burmester  constraint  pairs.  It  is  the 
author's  Intention  to  provide  a method  which  includes  the 
above  results  and  also  provides  the  capability  of  synthe- 
sizing a coupler  link  cr  ground  link  while  satisfying  five 
prescribed  multiply  separated  precision  points  for 
Burmester  constraints. 

A similar  algorithm  for  synthesizing  seven  multiply 
separated  precision  points  is  presented  for  the  generalized 
conic  pairs.  It  might  be  noted  that  for  conic  constraint 
pairs,  one  can  only  synthesize  the  coupler  link  since  conics 
(excepting  circular  constraints)  have  no  physical  center 
physically  represented  by  a fixed  pivot. 

One  might  say  that  the  dependence  of  dynamic  synthesis 
lies  solely  within  this  area  of  research  since  one  must  know 
at  least  one  constraint  pair  to  synthesize  the  absolute 
dynamic  characteristics.  This  particular  area  of  research 
represents  a distinct  discipline  of  kinematic  synthesis  and 
consequently  diversification  in  kinematic  utilization. 


kinematic’ an  more  versatile  in  his  ability  to  solve 
problems . 

This  paper  will  provide  an  algorithm  for  synthesizing 
seven  coplanar  positions  employing  conic  constraint  pairs. 
The  primary  emphasis  of  this  work  will  not  be  to  develop 
an  algorithm  for  synthesising  point-positions  using  conic 
constraints,  as  is  the  case  of  Burmester  contraints,  due 
to  the  complexity  of  such  an  endeavor.  The  author  will, 
however,  provide  avenues  to  such  research  and  therefore 
provide  a foundation  for  future  research. 

Thu  Application  of 

The  areas  of  application  of  coplanar  linkages  have 
recently  become  very  broad  and  divergent.  This  is  primarily 
a result  of  requirements  imposed  upon  design  by  competition 
and  society  to  create  safe,  reliable,  efficient  and  inex- 
pensive mechanical  systems.  This  has  brought  about 
considerable  demand  for  kinematicians  and  their  research 


Present  research  in  the  application  of  coplanar 


mechanisms  generally  falls  into  one  of  the  following 
categories: 


Multiply 


Generation 


( Positions 
J Trajectories 
\ Curvature 
I Stationary  Curvature 


Coplanar  ] (2)  Function 
Motion  | Generation 


Input-Output 
Relative  Displacement 


Multiply 

Separated 

Positions 

Coplanar 


Figures  (1-2... 1-7)  represent  a small  portion  o(  the  many 
applications  of  such  mechanisms.  These  applications  are 
self-evident  and  require  no  further  elucidation  at  this 
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INPUT/  OUTPUT 


Figure  (1-3)  Function  Generator  (Decimal 


hoy  Converter) 


CIRCULAR 

APPROXIMATIONS 


MATERIALS  HANDLING  BELT 
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^ * 1 With  Adjustable  Amplifier) 


CHAPTER 


BURMESTER  T 


Fundamentally  there  is  only  one  constraint  function 
used  by  the  kinematician  to  constrain  the  moving  plane. 

As  one  would  gather  from  Chapter  X,  this  is  the  conic 
constraint  function.  The  reason  is  that  for  convenience 
kinematicians  almost  always  employ  functions  which  are 
easily  constructed  physically  and  provide  utility.  These 
requirements  dictate  the  use  of  conic  constraint  functions. 
It  might  be  mentioned, however,  that  one  may  use  any 
analytic  function  as  a means  of  contraining  the  moving 

For  convenience  kinematicians  generally  separate 
circular  constraints  from  conic  constraints  since  the 
analytics  for  the  latter  are  significantly  more  complex. 
Figures  (2-1,  2)  represent  what  will  be  referred  to  as  the 
two  classes  of  mechanisms.  In  theory  the  two  classes  are 
the  same j however,  the  Burmester  Theory  represents  a 
degenerate  study  of  the  conic  constraint  algorithm. 

Notice  the  basic  difference  in  the  fixed  pivots  and  c 


*An  analytic  function  is  a function  that  can  be 
issed  by  a Taylor  Series. 


COUPLER 


PARABOLA 


(2-2) 


,))  ! ( (u2,  V,),  (H,,  Ka)}. 


EACH  POSITION  REQUIRES  3 PARAMETERS  {0,0 ,y) 


Figure  (2-3)  Variables  v-»-  Parameters  for 
Coplanar  Synthesis 


the  system  by  translating  the  pole  Poi  to  the  origin 
Figure  (2-4)  removes  parameters  ei  and  Bi-  rotating 
the  plane  and  scaling  so  that  e2  =*  1 and  G ; - 0 (Figure 
2-5),  leaves  the  following  parameters 


for  the  problem  specifications.  Hence  it  is  seen  that 


number  of  variables  = number  of  parameters . 

Therefore  it  is  shown  that  the  system  becomes  a functionally 
closed  system;  i.e.  a system  which  is  in  agreement  with 
the  variables  vs.  parameters  requirements. 

A simple  equation  governing  all  coplanar  constraint 
systems  can  be  expressed  as  follows 

dumber  of  positions  - number  of  variables  r 7 . 


Hence,  one  can  proceed  to  synthesize  Five  Multiply 
Separated  Positions  in  Coplanar  Motion. 


Using  the  cartesian  frame  of  reference  (rectangular 
coordinate  system) , the  transformation  of  the  coordinates 
of  ft(o,v)  as  a point  in  the  moving  plane  H,  to  the 
coordinates  of  the  point  A(U,V)  of  the  fixed  system  2 is 
given  by 

U - ucosy  - vsiny  + e 


ranslation  :'o  the  Pole  Reference 


v4 


Figure  (2-5) 


Rotation  and  Scaling  of  die 
Fixed  Reference  System 


From  equations  (2-1)  it  is  evident  that  the  origin  e of 
E(u,v)  has  the  coordinates,  e(a,3)  in  !'.r  as  seen  in  Figure 


The  intention  of  this  section  is  to  introduce  the 
seven  Am®  coefficients  which  are  sufficient  to  define  the 
motion  for  five  prescribed  positions  ( Z ) constrained  by 
Durr, ’.ester  circles.  The  results  can  be  in  terms  of  either 
crank  and  rocker  or  slider  crank  systems.  Writing  the 
general  form  of  the  circle  constraint  equation  G in  £ gives 


G (0, V)  = Q„(U* 


2Q5V  + Q, 


(2-2) 


For  a multiply  separated  position  differentiation  with 
respect  to  the  position  parameters  Y evaluated  for  a 
specified  Y.  gives 


* (2-3) 

a functions 


Writing  the  difference,  - G0,  gives  the  generalized  form 
for  the  circular  constraint  function  as 


Figure  (2-6)  The 


IU,V)  = Q {0* 


: + v*  - v,+ 

Q,(aui  - 20  ) + 0 (2Vt  - 2V  ) = 0. 


d^k{aoos7+  gsiny) 
^(-aSinyh  BoosY, 
cosy  - 1, 


Y =Y4. 


(2-6) 


VV  + W + Vu*  +vV 


after  substituting  equations  (2-1) 


respective 


U and  V terms.  By  rearranging  equation  (2-7)  in  the  two 
distinct  forms 


(Ao 


A,tU  + + [a»tJ  (uQ,  * 

t-An^u  + As.lQi  + [A.jjU  + A61)Q2  = 0 


[A,^  + A,jU  + A,^v]Q0  + lA,^!  (uQ2  - vQi) 

+1  Aj^u  + A.^IQ,  + [fts4v  + Aj^IQj  = 0 

provides  two  determinants 

|(Aoj  + Al^U  + Aj^v)  A|^  (-A^v  + As^)  (Afc^U  + Afj^)j  3 0 

i = 1,2, 3, 4 


and  (2-9) 

|^a°t  * aiJtU  + A?iV^  a,*t  ^a3fU  + a5f^  ^3ZV  * a6f^|  = ® 

both  of  which  must  be  identically  zero. 

The  two  determinants  appear  to  yi.eld  characteristic 
equations  that  are  of  third  order.  This  is,  however,  not  the 
case  since  the  determinant  provides  a nullity  of  the  third 
order  term.  Thus  these  determinants  provide  two  conics  of 
u and  v in  the  moving  plane  which  can  be  expressed  as 


?nv4  + (Vlsu  + Via)v 
?*lVa  + (722U  + 72,)v 


Vi«  = 0 


(2-10) 


Via  = lA,^ 

*”-K 

Via  = | At^ 


Vaa  » l&i, 

’»*  = i»>i 

Vas  “ [flij, 
Vaa  — | Ao  ^ 


A'l  S‘J 
A‘i  A‘J 

AsH  Astl 

A*i  A*J 

a»»  A‘*i 

a»*  r64 


A««  a*4 

A5i  A)*l 

A,t  A4j| 

A=i  A‘J 

Asj  A6(.| 


|Aat  Ajt  ASj  AaJ 
|Aoj^  Aa^-  R,,j  R6^| 

K A.,  A,,  A., | 
|a2j(  Aat  Ai4  A(t| 

]Ao^  Aa^  A*£ 

|rH  A'i  Alt  fl‘il 


The  intersections  of  these  conics  are  the  Burmester  points. 
Eliminating  v from  equation  (2-10)  gives  the  quortic 


Eiu4  + Eau1  + Eju*  + Eau  + Es  - 


in  u fcr  the  Burmester  points. 


(2-12) 

coefficients  are 


E,  = (c— h)  - f (c-li)  (a-f ) + h (a-f) 1 

Ea  = (h-c) I2(p-d)  + f(b-g)]  + (a-f) [f (p-d) 

+ g (h-c)  * p(a-f)  + 2h (b~g) ] 

* + 2 (h-c) (q-e)  + (b-g) [f (p-d)  + g(h-c) 

+ h (b-g ) + (a-f  J [f  (q-e)  * g(p-d)  + q(a-f) 
* 2p('o-g)] 


(d-p) 


E,  = (b-g)  [f (q-e)  + g(p-d)  + 2q(a-f) 
+(q-e) [2 (p-d)  + g(a-f)] 

Bs  = (e-g)a  + g(q-e)  (b-g)  + q(b-g)! 


a = 7 11  7, 
b = 7,i  7 i 
c = 7,i  7, 
d - 7,i  7 i 
e = 7,i  7i 


f = 7n  7,, 


p(b-g)  ] 


The  four  solutions  to  equation  ( 
removing  the  term  v'  from  equations  ( 
coordinate,  vn,  is  given  as 

f (p-d) unt 


(a-f )u_ 


T5~gy~ 


are  un.  By 


(2-15) 


The  location  of  the  centerpoints  in  the  fined  plane 


can  be  related  from 


£ + Al  a°n  + Rz  Bvn 
A,t  ♦ Aijun  " RUvn 


(2-17) 


(2-.13) 


» = (D<Fi  - DaPi)| 
n (B,F*  - EJP,  |n 


. Vis,  V,S(  Vi,, 


3 by  Figures  (2-7,3 


= 0,  i=  1,2, 3, 4, 


V2S,  Vs e = 0 
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Figure  (2-8)  Lincar-Poin t Synthesis 


' f), “■•>3  ■ '), 


Since  the  identity 

sinY*  + cosy?  - 1 (2-23) 

must  be  satisfied  one  may  write  the  requirement  for  a 
single  precision  point  as 


From  equation  (2-25)  one  can  write  the  positional  matrix 
for  positions  j - 0,1, 2/3, 4/5, 6, 7, 8 as 

8x6  V 0 


(2-2 


From  equations (2-26,27)  it  is  evident  that  there  are 
eight  unknowns 

{(Uj,  vt);  (Hi#  K.)}  t ( (u^ , vs);  0«a , Ka)  > . 

Thus  one  can  v/rite  the  eight  precision  point  equations  and 
theoretically  satisfy  his  solution  for  nine  finitely 
separated  precision  points. 


CHAPTER 


THE  STATE  OP  THE  ART 

Early  works  by  Baris  (10),  Beyer  111]  and  Hackmueller 

(12)  have  provided  the  kinematician  with  the  capability  of 
using  conic  constraints  for  prescribing  displacements  of 
the  moving  plane.  Later  works  by  Sandor  and  Preudens tein 

(13)  and  Woo  and  Freudenstein  [17]  provided  an  extension 

of  the  earlier  authors'  works  that  had  employed  infinitesimal 
displacements.  A recent  publication  by  Freudenstein, 

Bottema  and  Xoetsier  (8)  provides  the  kinematician  with  a 
concise  method  of  synthesising  six  finite  positions  of  the 
moving  plane.  This  paper  also  treats  all  degenerate 
cases  of  the  conics  for  finite  synthesis  through  six 
position  theory. 

The  intention  here  is  to  formulate  a generalized 
method  for  synthesizing  six  multiply  separated  positions  of 
the  moving  plane  by  defining  the  conic  section  point 
curve.  It  is  also  the  author's  objective  to  present  a 
new  algorithm  for  synthesising  seven  multiply  separated 
positions  of  the  moving  plane.  This  algorithm  produces 
two  quartics  in  u and  v which  relate  the  16  conic  points 
for  conic  constraints. 
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Figure  (3-1)  Typics 


Conic  Theory 
Variables  Vs .Parameters 


Before  synthesizing 
advantageous  to  analyze  t 
number  of  parameters  for 
Figure  (3-2)  that  there  e 
with  two  generalized  coni 


conic  constraint  pair,  it  is 
s number  of  variables  vs  the 
ic  system.  It  is  seen  from 
3 fourteen  variables  associated 


( *°i  'V^  ' (H,'  [j'ji1'  (0i>> 

{(us,v2),  (H2,  Ka> , j J|.iJ  ;e,  j|«  J ^ ] , (08)). 

After  normalizing  and  removing  a ,3  ,a  and  B by 
translation,  rotation  and  scaling,  as  outlined  for  the 
Burmester  constraints,  the  number  of  synthaslzable 
positions  becomes 

number  of  positions  = number  of  variables  + 7 
number  of  positions  = 14  + 7 
number  of  positions  - 7. 

The  maximum  number  of  positions  (coordinates  and  angles) 
which  can  be  synthesized  using  conic  constraints  becomes 
seven.  Hence  one  can  proceed  to  synthesize  Seven  Multiply 
Separated  Positions  in  Coplanar  Motion  for  conic  constraint 


VA 


EQUATIONS  OF  TRANSFORMATION 
(TRANSLATION  AND  ROTATION) 

U » (V-K)SINO  + (U  -H)COS0 

V = (V  -•  K)  COSO  - (U-H)SINO 


Figure  (3-2) 


ertie 


Position 


Theory 


The  equation  for  a generalised  conic  with  its  center 
at  the  origin  can  be  expressed  as 


where  Z , , , are  defined  for  the  conics  as  illustrated 
in  Table  (3-1).  If  the  center  of  the  conic  is  not  at  the 
origin,  the  inverse  transformation 

U =*  (V  - K)s£nO  + (U  - H)cosO 

(3-2) 

V = (V  - X)  cosO  - (U  - ir)sin0 
give3  the  new  coordinates  as  shown  in  figure  (3-2)  for  an 
ellipse.  Substituting  equations  (3-2)  into  (3-1)  gives 

Z ((V2  - 2VK  + K2)sin20  t 2 (OV  - UK  - VH  + KH)sin0cosO 
+ (U2  - 2011  t H 2 ) cos 2 G } + 

Z { (V2  - 2VK  + K2)cos20  - 2(0V  - VK  - VI!  t KH)sinOcos0 
+ (U2  - 2UH  + H2)sin20)  + 

7,  {(V  - K)cos0  - (0  - H)sin0)  + Z =0. 

J * (3-3) 

Rearranging  and  collecting  terms  gives 
U2{Z  cos20  + Z sin’e)  + UV{2  (Z  - 7.  Jsinocoso) 

+0{-2H(Z  cos2g  + Z sin20)  - 2K(Z  - Z )sin0cos0  - Z sin0} 
+V2{Z  sin  2 8 + 7.  COS*0)  (3-4) 

+V{-2K(Z  sin2©  t Z cos2©)  - 2H(Z  - Z )sin0cos0  + Z^coso) 
+ (H!  (Z  cos2Q  + Z sin29)  + K2(Z  sin2©  >■  Z cos20)  t 
2KH(Z  - Z )sir.0cos0  + 7 (Hsing  - :;coso)  + 7 ) = 0. 
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If  the  first  position's  origin  is  taken  to  be  coincident 
with  the  origin  of  the  fixed  plane,  then  a = 8 = y =0 
and  the  transformation  equations 

v ° + S "Y  + g (3-5) 

become 


This  allows  removal  of 

(H* {^008*0  t Z Sill* 0)  * K*(S  »ii»*8+  Z oos20)  + 

2KR (Z^  - Z^JsinOeosG  + 2^ (HsinG  --  K.cosG)  + 2 ^ } 

from  the  conic  constraint  equation  (3-4)  to  give 
(U2  - u*)(Z  co320  + Z sin20}  + 

(UV  - uv){3(Z  - Z )sinGcosG}+ 

(0  - u)  {-2H(Zicos20  + Z sin2G)  - 2K(2  - Z IsinGccsG  - Z sin0}t 
(V2  - v2){  Zisin20  + Z cos20)t 

(V  - v)  {-2K(Zisin20+  Z^cos2G)  - 2H(2^  - Z^JsinGccsO 
+ Z cosG } - 0 . 

* (3—6 ) 

Subscripting  the  transformation  equations  (3-5)  provides 

U3  = "°“*J  ' V8lnYj  + *j  (3.7) 

Vj  --  csinyj  + vcosYj  + R - . 

And  making  the  following  substitutions 


= 2(Z>  - Z^)  sinOcosO 
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Substituting  the  transformation  equations  (3-7)  into 
the  first  typical  equation  of  (3-10)  will  give 

(8i3(u2  - v2)  - 23^  j (uv)  + 2Bsj(u)  - iB.jfv)  + Oa.j)  X(  * 
{B^fu*  - v2)  + ZB^Cuv)  + (B^j  + B^)  (u)  +(Bsj  - B^j) 

(V)  + Bitj}*t+ 

(Bj:j(u)  - B (v)  + + 

- v!)  + 2Baj(uv)  t 2B(j(u)  + 2B7j(v)  + BijS}X(  * 
{B^ j (u)  + B,j(v)  + B,,jJXs  = 0 

(3-11) 

For  any  multiply  separated  position  a typical 
equation  can  be  written  as  above  by  taking  the  coefficients 
Bm£  from  Table  (3-2) . For  a finite  position  the  coefficients 


0 S 

i-  f i:  v 

<# 

'i  s 1 1 !* 

a? 

p 

a? 

-fl  i rij? 

1 

a? 

H'H  ' 1 

a? 

flh-i  1 i 

a? 

•is"!  ! 1 

a? 

riHU  ! 

a? 

Si? 

r f f ■ • 

a? 

Iff  • • 

a? 

fit  ■■  • 

ara  given  by 


B.j  = "3inlTj 
B^  = siny -eosYj 
Bjj  = cosYj  - 1 


which  may  be  expressed  similarly  to  the  coefficients  of 
equations  (2-6) . Prom  the  typical  equation  one  can 
develop  the  Six  Position  Conic  Constraint  Matrix  for  six 
multiply  separated  positions  as  defined  by  positions 
i = 0,1, 2, 3, 4,5  of  the  moving  plane. 

The  Six  Position  Conic  Section 


As  shown  by  Freudenstein,  Bottema  and  Koetsier  the 
Conic  Section  Point  Curve  is  of  seventh  degree.  This  is 
verified  from  equations  (3-11)  by  adding  column  1 to 
column  4,  which  gives  for  the  typical  equation  after 
subscripting  for  ith  position 


{3zi<u2  - v2)  + 2Bit(uv)  + 

(2,ButV1(u|  +2(B7l- 

{Buw  + + b.»£}xs 


Vv 


>t  *£' 


(3  x 3)  in 


g^Uu*  - v*I*  * 2(»*  - v'luv  + 4(uv)*}  + 

<92mu  + 9jmv)  lu,_  v*5  + 

<9„n,“  + 9smv  + g(U>  {uv}  + (3-16) 

V*"*1  + + 

(V"  * 910rav  + »„■> 


V {4lB.p  B.»"  B.»V  - 2[V(B,p  ' V 

* W“  v» 

9S™  - t2[B,iV  - B,»B,.P  + W - V 


-B.i  (B«p + V1 

-21VB,1  - V,.*  + Bsp(B,i  - V 

-V‘v  vn 


V {tBuB,.P  + B,iV  -(B.pb..*  + W 
+ 21b,*(b«p  + V -B5P(B,t  + B^]1 
“ t<8»**UP  + B,*V  - (B.PB..l  + B,pB,l> 

+2,Be,(BSp  - V - W - B.*>1} 

V={2‘VB„p-BspB..t> 

+ B,*fBsP  + V>  " W V> 

g1.™“{2(B,pB1„-B.1B,.P> 


+B,*(B,P  - V>  - B,P(BS* 

{B,»B..p  ' VB..i} 


n = 1,2,3. ..10. 


sVl)3  Om[gim{(u!  - v*)*  + 2 (u!  - v*)uv  + 4 (uv) z) 

+‘VU  + Vv){u’  ' vM 

+ (9,»u  * 9S™V  + Vtuv} 

+ (VU  + 9,.mv  + *„»>!  =°- 


a + 2 + 3,  + it 


property. 


Appendix  [A],  the  cobles  Q„ 


(Himu  + Hjm)  (u2  + vz)  + Hjrau2  + 
H,mv!+  Hsmuv  + H,raU  + H,mv  + «„.«  " Om- 


(3-19) 


By  rearranging  the  first  term  of  equation  (3-18)  the 
highest  order  of  the  conic  section  point  curve  becomes 


2uv((uz  - v2)  + uvj) 

where  1 and  1 are  scalers  dependent 
8U1  ^ . it  can  be  shown  by  substituting 


(3-20) 

coefficients 


into  equation  (3-18)  and  letti.ng  vw»  that  the  resulting 
slope  M is  related  by 

(M2  + « * (X  M +X  ) = 0 

1 2 (3-21) 

which  will  give  six  imaginary  asymtotss  indicating 
tricircularity . 

The  Describing  Parameters 


For  a particular  u and  v (-.he  conic  constraint  matrix 
le]  becomes  defined.  The  application  of  Cramer's  Rule 
provides  for  a particular  u and  v 


Substituting  equations  (3-; 


(3-22c)  into  (3-8) 


EtL^cos8©  - 2sin©cos©)  (L  sin2©  + 2sin0cos©)| 
[L^cos8©  - sin8©)  (t  cos8©  - sin8©) 


By  adding  columns  1 t 
system  reduces  to” 


d dividing  1 by  cos8©  t 


1?.  - 


Substracting  column  2 from  column  1 yields  the  following 


( - 2Tan0)  (L  ) I 

[(1  - Tan8©)  (^-1)1 


2 determinant, 


a angle  of  rotation  function. 


The  centerpoint  coordinates  are  given  by  equations 
(3~8c)  and  (3-3e)  for  Z /z  = 0 ss  taken  from  Table 


(3-1)  for  the  circle,  ellipse  and  hyperbola 


Substituting  of  equations  (3-22)  gives 


(3-26) 

which  allows  for  expressing  the  centerpoint  coordinates 


becomes 
a,{(v‘  -: 

V<vz  " 


letting 
Y|  - (Vs 

Y,  = (yS 


! = 0,  the  constraint  equation  (3-3) 

■ K!)  sin!Q  + 2(uv  - UK  - v!I  t KH)sinOcos0 
(ul  - 2uH  + az)cosz0l  + 

+ K1 ) cos!Q  - 2(uv  - uK  - vH  + KH)sin0cos0 
(u*  - 2uH  + Ha)sin*0)  + Z = 0 


; + Kz)cos*0  - 2(uv-  uK  - vH  + KH)sin0cos0 
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and  combining  with  equation  (3-23b)  gives  the  system 


The  inversion  provides  the  semi-axis 

Zi  - -(Ljsin'O  - cos -Q)  

Z,,  VifhjSinaG  - cos *3)  - yt(L7oos>e  - sin2©) 

Z 5 = (L,cos*3  - sln2Q) 

z*  Y, (L3sinJ0  - cos'S)  - (lscos!©  - sin'O) 

(3-29) 

All  of  the  properties  for  conic  constraint  synthesis  of 
six  multiply  separated  positions  can  now  be  expressed  for 
the  circle,  ellipse  and  hyperbola.  As  will  be  shown  in 
the  following  section,  the  parabola  requires  a different 
algorithm  in  locating  the  center  and  angle  of  rotation  for 
the  semi-axis.  The  geometric  center  of  the  parabola  is 
at  infinity  which  requires  that 


A method  of  relating  the  center  as  the  vertex  is  derived  and 
illustrated  for  this  special  conic. 

The  Describing  Parameters 


For  a particular  u and  v the  conic  constraint  matrix 
[C]  becomes  defined  and  the  application  of  Cramer's  Rule 
provides 


for  the  particular  set  of  coordinates  (u,v) . From 
equation  (3-1)  and  Table  (3-1)  if  the  condition 

I.  * — «.  « 0 (3-31) 

is  satisfied  then  the  conic  is  a parabola  and  Table  (3-1) 
suggests  that 


For  a parabola  equations  (3-8)  become 


X = ZZ^inGcosG 
X = -2HX  - KX  - Z sinS 
X = Z sina0 


s gives  for  equations  (3-3 


(3-32) 


both  satisfying  equation  (3-31) . The  angle  of  the 
transformations  can  be  expressed  as 


It  will  be  shown  later  in  this  text  that  determines 

the  direction  which  the  parabola  opens  and  establishes  the 
phasing  of  0.  Equations  (3-32c)  and  (3-32e)  can  be 
grouped  to  form  the  characteristic  set 


The  discriminate  of  equation  (3-34)  can  be  expressed 
by 


as  defined  by  equation  (3-31) . If  equations  (3-34)  are  to 
relate  the  coordinates  of  the  center,  they  must  be 
coincident.  The  coincidence  also  requires  the  discriminate 


to  be  zero.  The  discriminate  gives  for  z /Z 
Z,  = cosed. iLt  - 2L3L;)  . 


Li  (1  - L,) 


B{  {-2v*i  + uXa>  + 


Bn  {X  } + 

B,;  OO  + 

Bw  {:<S}  = ° 
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(3-42) 


which  comes  iron  equation  (3-11) . For  seven  positions  of 
the  moving  plane  the  system  can  be  expressed,  by  applying 


where  o^  are  functions  of  the  Bm(_  coefficients. 

The  first  two  expressions  of  equation  (3-43)  are 
functions  of  only  X , X , X and  may  be  expressed  as 


((u*  - v4)  - 2e(iu  + e v - e^JX  + 

{ 2uv  - (e]  + e lu  - (e  - e^)  - e )X  + 


{ -(u4 


Uus  - v1)  - (e  + e s)u  - (e  - e )v  - e t>X 
(2uv  - e u + esjv  - e^JX^  = 0- 


(3-45) 


Though  these  two  equations  were  linearly  derived  from  t 
typical  or  generalized  set,  they  must  satisfy  the 
characteristic  set  and  the  characteristic  equation.  By 
multiplying  equation  (3-44)  by  Bj  and  subtracting  from 
the  2th  equation  (3-11)  and  by  multiplying  equation 
(3-45)  by  , end  subtracting  from  the  2th  equation 
(3-11)  gives  for  the  2th  or  typical  equation. 

(f  u t f v + f )*  + 


columns  4 and  5, 


coefficients  of 


are  not  altered  by  this  manipulation.  Employing  the 
technique  outlined  for  the  Five  Position  Burmester 
Theory  for  compound  determinants  with  linear  orthogonal  or 
linear  parallel  elements  gives 


(3-49) 


The  two  determinants  (characteristic  equations)  of  the 
square  matrices  provide  two  quartics  in  u and  v relating 
the  16  Burmester  points  for  the 
can  be  expressed  as 


general 


quartics 


where  the  4pq's  are  the  expansions  of  the  respective 
determinants . The  ipq's  can  be  expressed  as 


and  similarly  for  equation  (3-51) . Thus  the  quartics 
defining  the  16  Burmester  points  can  be  related  by  a closed 
form  technique. 

The  intersections  of  these  quartics  can  be  found  by 


6.1 


“ *n»  + Am<v  + Ara>v  + AraioV 
^m5  = 4mi  i + AmiaV  + AmnV  + 


and  employing  Sylvester's  Dialytio  Method  of  elimination 
to  give  the  resultant  R for  vn  as 


For  the  real  Burmester  points  the  resultant  R(vn)  must  be 
singular.  Therefore  various  numerical  methods  of 
minimizing  R may  be  employed  to  obtain  the  real  Burmester 
points  (maximum  of  16) . 

The  describing  parameters  of  the  conio  may  be  found 
using  the  same  technique  outlined  for  six  position 
synthesis  equations  (3-22, .. .3-41)  for  each  of  the 
intersections  of  the 


I' tics  (Burmester  points) . 


CHAPTER 


ANALYSIS  OF  THE  GENERALIZED 
CONIC  CONSTRAINT  PAIR 

As  shown  by  Sandor  and  Freudenstein  [13] , synthesis  is 
not  and  should  not  be  restricted  to  circular  (Burmester) 
constraints.  In  general,  any  conic  constraint  set 
provides  a means  of  constraining  the  moving  plane  which  is 
the  requirement  of  synthesis.  This  is  to  say,  and  will 
be  verified  in  what  follows,  that  the  Burmester  constraints 
are  merely  degenerate  conic  constraints  which  utilise 
simplified  modes  of  synthesis  and  analysis.  Since  the 
researchers  Sandor  and  Freudenstein  [13],  Freudenstein, 
Bottema,  and  Koetsier  [8]  and  this  author  have  algorithms 
for  synthesising  positions  in  coplanar  motion  for  conic 
constraints,  it  seems  only  logical  that  an  algorithm  of 
analysis  is  necessary.  Therefore  it  is  the  author's 
intent  to  introduce  such  an  analysis  algorithm  for  any 
combination  of  the  following  conics 

Circle 

Ellipse 

Hyperbola 

Parabola. 

From  equation  (3-6)  and  equation  (3-8)  the  conic  constraint 
equation  was  given  as 

(0*  - U!)X  + (BV  - uv)X  * (U  - U)X  + 


(V  - 


1-1) 


= ucosy  - vsiny  + a 


= 2(Z  - Z JsinOcosO 


z sin*©  + Z eos’o 
-HX  - 2KX  + Z COS0. 


a known  constraint  s 


become  known.*  Rearranging  gives  the  quadratic 
A a8  t(A  I t 1 lo  t 1J!  * M * J = 0 


! siny{X  u - 2X  v}  + cosy{2X  u + X v)  + 

= siny{2X  u - X v}  + cosy { X u + 2X  v)  + 
• -sin8y{ (u8  - v!) (X  - X ) + 2uvX  ) + 
sinycosy{(us  - v8)X  - 2uv(X  - X )) 


(cosy  - 1) (uX  + vX  } + 


equations  (4-5)  nay  be  subscripted  as  fellows 

(xm,n) 

n ^ 1*2, 3,4, 5, 6, 


The  two  constraint  sets  prov5.de  two  quadratics  for 
equation  (4-4)  and  may  be  expressed  as 


(4-6) 


in  terms  of  y for  the  moving  plane.  Applying  Cramer's 
Rule  to  equation  (4-6)  squaring  the  second  row  and  equating 
to  the  first  row  gives  the  quartic  for  6 


(4-7 


HYPERBOLA 


in  terms  of  y of  the  moving  plane.  The  roots  of  the 
quartic  provide  four  8 coordinates  of  the  moving  for  each 
angle  y provided.  The  corresponding  a is  taken  from  the 
Inversion  of  equation  (4-6) . The  coefficients  E^  are 
defined  by  equations  (2-13)  and  (2-14). 

Higher  Order  Contacts 

It  is  often  advantageous  to  analyze  the  geometric 
velocity,  acceleration,  jerk,  jounce,  etc.,  of  a particular 
constraint  pair.  This  is  accomplished  by  differentiating 
the  Position  Matrices  equations  (4-6)  with  respect  to 


Geometric 

Velocity 


fm,  “ 2Am,°  * Am2B  + 4ms 


(4-9) 


Differentiating  of 
acceleration 
Geometric  f 
Acceleration 


Velocity  Matrices  will  give  the 


Matrices 


(-10) 


Differentiating  equation  (4-10)  provides  the  jerk 


This  illustrates  that  differentiating  and  applying  Cramer's 
Rule  provide  a mode  of  analyzing  all  orders  of  contacts. 


e higher  order  contacts  are  dependent 


lower  order  contacts. 


(A  -A  ) 8 +(A 


CURVE 


Coupler  Curve  Function  for  Burmester 
constraints 


greater  versatility  in  synthesis  using  conic  constraints 
as  has  been  proven  to  be  the  case. 

Locating  the  Pole  of 


The  pole  or  instant  center  is  the  position  in  the 
moving  plane  which  has  an  instantaneous  velocity  of  zero. 
And  since  the  pole  coordinates  are  given  by  the  trans- 
formations 


differentiating  with  respect  to  y gives  for  Up  = Vp  = 0 


= Upsiny  + vpcosy 
S ' = -Upcosy  + VpSiny. 


(4-13) 


Substituting  equations  (4-13)  into  equations  (4-8) 
provides 


Applying  Cramer's  Rule  gives  the  pole's  position  in  the 
fixed  plane  for  any  y chosen.  Prom  the  pole  position  the 
location  of  the  inflection  circle,  Ball  point  and  other 
can  be  determined. 


higher  order  properties 
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(4-17) 


hohe  that  the  sign  of  Z /Z  is  determined  by  G and 
establishes  the  phasing.  For  the  initial  position  '-'  0 , 
the  constraint  equation  (3-4)  becomes 


H!(Zicos!0)  + 

KM^sin'G)  + 

KU(22  sinGcosO)  + 

H(Z  t-inG)  + (3-37) 

(uaX  -V  uvx  + u*  t v’x^  + vX  } = 0. 

Multiplying  equation  (3-32c)  by  H and  multiplying  equation 
(3-32e)  by  X and  summing  the  two  relations  gives 


K!  {■/,  coszG)  + 
KJ(Z!S.in*0)  + 
KB(2Z  sinGcosO)  + 


K (Xs  - ZiCOSp)  _ 0_ 

2 (3-38) 


Sul>traeting  the  above  equation  from  equation  (3-38) 
provides  the  linear  expression 

(X  - 2 sinGllI  + (X  +7.  cosOlK  = -2X 

where  (3-39) 

X = {uzX  + uvX  + uX  + v*X  + vX  1. 

Combining  this  equation  with  (3-' 


The  inversion  of  the  equations  provides 

H = -4X,X6  - ( -X 5 * Z;.cos0)  (X3  + ZjCOsO) 

2X«i  (Xs  - ZssinO)  -Xi  (Xs  + Zjsin9) 


K = (X,  - Z,sin0) (~XS  + ZjCosO)  4 

Simplifying  gives 


[I.J 'sino  t-  2cosQ]  j'zj'p'  HsinioX,  - 2cos'ol,l| 


X = - L [z|]  + 2 [L^cosO  * J.us 


- 2 [t  LiicosJ0-L  [z6J 


(L  “sin©  2cosG] fzJ1  - (4sin*0L  - 2cos*0L  ] 


e above  coordinates  relate  the  position  of  t 
r the  specific  parabola  corresponding  to  the 
d complete  the  description  of  the  parabola  i 


Seven  Position  Conic  Theory 


In  the  previous  section  it  was  shown  that  six  positions 
of  the  moving  plane  provide  a seventh  degree  polynomial 
in  terms  of  u and  v (coordinates  of  the  conic  points)  in 
the  moving  plane.  The  addition  of  a seventh  position  will 
provide  a second  seventh  degree  polynomial  in  u and  v 
indicating  that  there  are  49  intersections  of  the  two 
curves.  As  shown  by  Freudens tein , Bottema  and  Koetsier, 
of  these  49  theoretical  solutions  one  must  subtract  18 
for  the  circular  points  at  infinity,  10  for  the  common 
poles  and  5 for  the  common  Ball  points, leaving  16  Bur- 
mester  points  for  the  general  case.* 

To  obtain  the  determinant  of  the  5x5  matrix  in 
terms  of  u and  v is  a formidable  task,  therefore  it 
becomes  desirable  to  express  the  matrix  in  some  other 
format.  If  the  system  is  arranged  in  terms  of  common 
coefficients,  the  generalised  equation  of  motion  (typicial 
equation)  can  be  expressed  as 

B ((u*  - v*) (X  - X ) + 2uv(X  )>  + 

B ( (u*  - Vs)  (X  ) - 2uv(X  - X )}  + 

Bj  {uX  + VXs)  + 

B>  (-VX  + U*s)  + 

Bs  (2uX[  + vX^)  + 


*The  circlepoints  at  infinity  are  a result  of  the 
function  being  tri-circular  as  previously  shown  in  this  text. 


velocity,  acceleration,  jerk,  etc.,  of  every  point  in 
the  moving  plane. 


Conversion  to  Seal  Time 

Emphasis  here  should  be  made  that  conversion  to 
real  time  requires  that  differentiation  be  taken  with 
respect  to  time  and  not  the  angle  Y • The  absolute 
velocity  can  be  expressed  as 

Absolute  I do  = X dy  + do 
Time  I dt  ‘dt  dt 

Velocity  < (4-19) 


For  the  absolute  acceleration  one  simply  differentiates 
equation  (4-19)  to  obtain 


Absolute 

Acceleration 


dte  = X i Hta  - X 


and  similarly  for  higher  order  motions. 
Observing  equations  (4-19,20)  it  is  s 


relations 


time.  By  judiciously 


assuming  the  origin  at  the  circlepoint  o£  constant  or 
defined  id  (rate  of  angular  rotation)  then 


R =*  centerpoint  to  circlepoint  length 
Rp  = circlepoint  to  pole  length. 

For  higher  order  derivations  differentiation  must  be 
taken  with  respect  to  time  for  a,  B,  Y.  u,  Rc  and  I?p. 


CHAPTER 


LINKAGE  SYNTHESIS  FOR  THE 
GENERALIZED  CONIC  CONSTRAINT  SET 


Often  the  kinematician  must  have  the  ability  to 
synthesize  a conic  constraint  set  while  satisfying 
coordinate  positions  (precision  points)  in  coplanar 
motion.  This  requires  forfeiting  the  angle  y of  the 
moving  plane  or  the  absolute  values  of  the  derivatives. 

From  Chapter  III  the  generalized  conic  constraint 
equation  was  given  by 


+ (V!  - vz)X4  t 2 (V  - v)Xs  =■  0 


(5-1) 


where  U and  V are  defined  by  equations  (3-7)  and  X ■ by 
equations  (3-8) . specifying  the  following  parameters 


and  letting  each  position  carry  the  subscript  i will 
give  for  the  finite  case 


siny^  t siny^cosy^  * a^  siny^  + a^cosy^  + a^s  - 0 


r r r” 


a.,  = (-(Xi  - V(u>  - v’)  - 2uvX,} 

a.a  = {Xa(u!  - v2)  - 2<X(  - X^uv) 

ais  = - 2vXj)  - Si(vXi  - 2uXt)  + (uXs  - vXj)} 

ai,  " <“i<2uX,  + vV  + 8i(uX,  + 2vV  + (UX3  + VX*>J 

ai5  - foiXi  + Hhx,  + aix,  * ®ix%  + 8ixs  - 'uX3  + vXs>>- 


citsin‘Yi 


ciisin>''i  + Ci,sin2Yi  + e^sinYi 


Cij  = 0 

(5-4) 


= ai!  + ai. 

= 2 (aj.1ax1  + •i,*!,) 

3 ai,  + ai l ~ ai,  + Zaitals 
“ 2(8i3ais  ' ai,ai,J 

= ail  - av 


(5-5) 


a*  and  S* 
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=.  a 

r,  r2  % 

0 

0.0  0.0 

0®  0®  0®  0® 

« 

.3  -.6 
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Conic  Set  Synthesis  While  Satisfying 
Seven  Precision  Points 


where  the  f^'s  are  taken  from  Table  (4-1).  In  synthesizing 
infinitesimal  positions  the  kinematician  can  specify  only 
the  relative  change  of  the  invariants  by 


Expressing  equations  (5-6)  and  (5-7)  as  a linear  set  for 
tile  first  derivative  gives 


(f^Jja*  + (f  M'  = -<f,,>i 


Applying  Cramer's  Rule  provides  the  proper  a' , 6'  for  the 
order  contact  synthesized.* 

For  higher  order  synthesis  one  need  only  differentiate 
equation  (5-8)  further  as  given  by  equations  (4-10,  11,  12) 
take  the  coefficients  fpg  from  Table  (4-1)  and  proceed  as 
outlined . 


Linkage  Synthesis  for  the 
Burmester  constraint  i'air 

As  one  might  gather  from  the  previous  chapters, 
Burmester  pair  synthesis  is  an  elementary  study  of  conic 
section  constraint  synthesis.  From  equation  (2-5)  the 
Burmester  constraint  equation  is  given  by 

(U1  - u2  - V2  - vs)Qt  + 2(U  - u)Qi  + 2(V  - v)Qt  = 0 


*These  are  the  unique  absolute  values  to  be  used  in  the 
synthesis  algorithm. 


where  U and  V are  taken  from  equation  (2-1)  and  Q /Q  from 
equation  (2-16) . By  specifying  the  following  parameters 

Constraint  Pair  I 

and  { (u,  v,  H,  K,  a,  8) 

Precision  Point  I 

and  letting  each  position  carry  the  subscript  (i)  gives 
for  the  finite  case 

aiisinYl  + a^cosy.  + atj  = 0 (s_9) 


which  is  a transcendental  equation  of  second  degree  in 
y.  The  coefficients  are  expressed  by 


“ (a,  - H)u  + (8.  - K)v 


Rearranging  and  squaring  provides  the  quadratic 

°i1stn*y  + ci,alnT  + ci,  - 0 


i-12) 


Figure  (S-2)  Burmestor  I.ink  Synthesis  While  Satisfying 
Five  Precision  Points 


also  exists  for  the 


The  uniqueness  of  a and 
Burmester  constraints  for  the  relative  change  prescribed. 
Differentiating  equation  (5-9)  gives  for  the  first  order 
infinitesimal  case 


where  the  f^  coefficients  are  taken  from  Table  (4-1) . 

In  synthesizing  infinitesimal  positions  we  specify 
only  the  relative  change  of  the  invariants 


Combining  equations  (5-13)  and  (5-14)  gives  the  linear 
system 


(f^JjV  + (f12>iB'  = 


Applying  Cramer's  Rule  provides  the  proper  a'  and  3'  for  the 
order  contact  synthesized.* 

For  higher  order  derivatives  one  need  only  differentiate 
equation  (5-.9)  further  as  given  by  equations  (4-10,  11, 

13) , take  the  coefficients  from  Table  (4-1)  and  proceed  as 
outlined. 


e corresponding  f q 


Coupler  and  Fixed  Link 
Synthesis  for  Conic  constraints 


S3 


It  often  becomes  desirable  to  synthesize  the  pin 
coordinates  of  the  coupler  link  or  fixed  pivots  of  a 
conic  constraint  pair  while  satisfying  seven  multiply 
separated  precision  points  in  coplanar  motion.  The 
designer  is  often  faced  with  such  requirements  when  the 
solution  to  a problem  dictates  size  limitations  or  optimal 
constraint  locations. 

It  is  important  here  that  one  investigate  just  what 
can  be  prescribed  while  satisfying  seven  precision  points 
in  coplanar  motion.  For  seven  precision  points  there  are 
six  angles  y.  which  indicate  that  there  are  six  prescrib- 
able variables  for  the  seven  precision  point  study.  For 
higher  order  studies,  the  angles  would  be  replaced  by 
the  relative  change  of  the  invariants 


as  explained  earlier  in  this  chapter  and  given  by  equation 
(5-7) . The  six  prescribable  variables  may  be  any 
combination  or  variation  of  the  following 


{(u.v^, 
((u^v^, 
{(ui,vi), 
{(K  ,H  ), 


(0^0  ),  (0,8)} 


(0t,  0 ),(a,B)). 


(a,  8) ) 


Thus  one  has  gained  by  the  versatility  o£  specifying 
constraint  requirements  while  satisfying  seven  multiply 
separated  precision  points. 

The  generalized  conic  constraint  pair  equation  (5-'!) 
is  rewritten  here  as 

ct  sin'Yi  + sin’y^  + c^sin'yj  + ei(sinyi  + cis  = 0 

(5-17) 


as  defined  by  equations  (5-5) . Unlike  the  conic  constraint 
pair  synthesis,  the  study  herein  is  concerned  with 
synthesizing  variables  for  both  constraint  pairs.  Therefore 
for  each  precision  point  i there  exist  two  functions  of 


Constraint  #1  (f^sin^  + £ijSin3yi  + fijsin!yi  + ^sihYt 


Constraint  82  {gj^sin'yj  + gj^sin’y^  + g£jstn!y^  + gj^siny^ 


as  taken  from  equations  (S— 4)  and  (5-5) . 


resolution  of 


fcrigancnetric  relations 


seems  to  be  a formidable  task.  And  since  the  algorithm 
requires  the  sets  of  fp(J  and  gog  which  have  a common 
root  or  roots,  a mode  of  solution  becomes  necessary. 
Employing  Sylvester's  Dialy tie  Method  gives  the  resultant 


3i 


3iS 

(5-19) 


This  is  to  say  that  the  determinant  R^  must  be  singular 
for  each  of  the  precision  points  prescribed.  Therefore 
one  can  prescribe  any  six  of  the  following  variables  in 
addition  to  the  coordinates  of  the  seven  precision  points 


(uj , a,,  kj( 

lV  V V V 


(5-20) 


by  employing  a numerical  search  routine  satisfying  equations 
(5-19). 

For  the  infinitesimal  positions  equations  (5-16) 
and  the  derivatives  of  the  constraint  equations  (5-2) 


Constraint  {1  {lf’V)1aIt  + + tl'^) j_  = 0 

Constraint  S2  {(f^)^a!c  + lf^)j8*  + (f^)^  = 0 

Prescribed  { b.  o'4  + =0 

Orders  1 8 


(5-21) 


where  the  coefficients  fmn  are  taken  from  Table  4-1. 

Therefore  the  requirements  for  any  order  of  K 
become 


Thus  for  any  finite  or  infinitesimal  precision  point  i, 
the  requirements  imposed  by  the  prescribed  precision  point's 
coordinates  and  the  conic's  prescribed  variables  equations 
(5-20)  are  given  by  equation  (5-19)  or  (5-22)  depending 
on  the  case  study. 


Coupler  and  Fixed  Link 
Synthesis  for  Burmester  Constraints 

As  one  may  surmise  from  previous  discussion,  the 
synthesis  of  coupler  and  fixed  links  for  the  Burmester 
constraints  represents  a simplified  form  of  the  conic 
constraint  algorithm.  Rewriting  the  generalized 
Burmester  constraint  equation,  equation  (5-9) , for  the 
finite  case 


where  a^j  are  defined  by  equations  (5-10).  Prescribing 
variables  for  both  constraints,  results  in  the  two 


Constraint  SI  {gijsinYi  + gtjcosYi  + gis  = 0 
Constraint  S2  {fj^sinYj  + cosYj  + ftj  = 0 


gij  = f*u,'  v,<  H,' 

fij  = f(uj(  vt,  Hj(  Kj(  a,B> 

i = 1,2, 3, 4 
j = 1,2,3 


as  taken  from  equations  (5-10).  Equations  (5-24)  require 


By  synthesizing  the  coupler  link  variables 
{ (u  , v ) , (u  , v ) } 


will  give  for  the  unknown  variables  of  equation  (5-26) 


» U-. 


aij  = [B1 


1x6 
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By  specifying  the  fixed  link  variables 
{(Ht,  K(),  (Hjf  Kt)) 
provides  from  equation  (5-26) 

rci 

5x6  u v 


(5-29) 


where 


Therefore  one  need  only  satisfy  the  above  requirements 
for  synthesis  of  either  the  coupler  link  or  fixed  link. 

As  for  the  infinitesimal  displacements  the  procedure  is 
identical  to  that  outlined  by  equations  (5-21)  and  (5-22) . 


CHAPTER 


SIX  MULTIPLY  SEPARATED  POINT-POSITIONS 
IN  COPLAHAR  MOTION 

In  Chapter  II  the  generalized  Burmester  equation 
(2-7)  was  given  as 

A Q * Q^u  + Aj  Qv  + Aj  (uf^  + vQ^) 

+ A (uQ  - vQ  ) + A Q + A Q = 0. 

”t  » 1 5i  1 H 1 

Making  the  following  substitutions 

Q = Z (uQi  + vQa)  - 

Q u = Zj  (uQ2  - vQt)  = Z^ 

Q v = Z Ql  - Zs 

gives  the  linearized  expressions 

! A Z„ , = 0 
m=.0  mZ  “ 

(6 

as  the  required  circular  constraint  formulation.  The 
six  point-positions  of  the  moving  plane  are  Z - 0.1,2, 
where  the  first  position,  l - 0,  is  superimposed  h 
fixed  origin.  The  last  five  of  these  positions  may  be 
expressed  as 


[A]  (2)  - 
5x7 


(6-4 


In  prescribing  six  point-positions  of  the  moving 
plane,  as  illustrated  in  Figure  (6-1) , the  angular 
rotation  of  two  of  the  positions  must  be  compatible 
so  that  all  combinations  of  five  positions  have  two 
common  Surmestar  pairs.  Similarly,  infinitesimal  cases 
require  invariant  magnitude  compatibility. 

Theorem  I 

There  exists  a unique  pair  of  Burmester  constraints 
(pair)  for  six  multiply  separated  position  in  eoplanar 
rotion . 

For  six  multiply  separated  positions,  equations 
(6-4)  can  be  reconstructed  in  terms  of  3 unknowns.  This 
can  be  accomplished  by  Gauss-Siedel  reduction  or  by 
Cramer's  Hula.  Tills  allows  one  to  establish  the  two 


•J,  h 
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Figure  (6-1)  Six  Multiply  Separated  Point-Positions 
in  coplanar  Motion 
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The  linear  reduction  requires  that  all  of  the  Bu raves ter 
points  must  satisfy  the  line  defined  by  equation  (6-5a) 
and  all  of  the  Burmester  centers  must  satisfy  equation 
(6-5b) . from  Chapter  XI  it  was  shown  that  five  positions 
uniquely  determine  the  Burmester  points  as  the  inter- 
section of  two  conics  (the  centers  may  also  be  determined 
by  two  conics) . Since  the  lines  defined  by  equations 
(6-5)  can  only  intersect  a conic  twice,  there  exists  a 
unique  pair  of  Burmester  constraints  satisfying  six 
multiply  separated  positions  in  coplanar  motion. 

Algorithm  for  Six  Multiply 
Sepa rated  Point-Positions  in  coplanar  Motion 

For  simplicity  the  following  algorithm  will  be 
presented  in  tha  most  general  sense,  i.e.  it  i3  given  in 
terms  of  the  requirements  for  all  cases.  The  formulation 
will  show  that  each  case  study  requires  a different  set 
of  analytical  expressions  to  define  the  compatibility 
angles  or  the  derivative's  magnitude  depending  on  the 
case  being  studied. 

Rearranging  equation  (6-3)  into  the  form 


(6-6) 


3nd  applying  Cramer' v 


jives,  after  a sign  adjustment, 


3-9)  gives  bi- 


quadratics 


(6-10) 


{V(?s] 


(?.2VS1 


(V  7 


e!l  “ "(VV3I  + ?,,V,  1 • 


+ 7M7Jt)} 


+ 7 V )} 


(6-11) 


Prom  equations  (6-2),  the  identity 


shows  that  equations  (6-10)  can  be  written  as 


Equations  (6-12)  must  have  the  same  roots  as  determined 
by  Theorem  I so  that 

IB]*  3 0 


(6-13) 


®pq  = f(fl_  ) 

>■  ra  = 0,1,2, 3, 4, 5, 6 

1 = 1,2, 3, 4, 5. 

This  allows  one  to  express  two  singularity  constraint 
functions  as 


(6-15) 

••J 


depending  on  the  case  of  study. 

Once  these  singularity  constraint  functions  are 
satisfied  the  quadratics  equations  (6-12)  become  defined 
and  the  Burmaster  points  are  related  by 


given  by 


and  the  v coordinates  a 


s taken  from  equations  (6-7) . The  coordinates  of  t 
enters  are  related  by 


“n  = -K, 


Thus  one  has  gained  the  versatility  of  specifying  four 
ooplanar  positions  and  two  precision  points  for  the 
moving  plane.  This  study  will  be  refered  to  as  a point- 
position  study  since  it  is  a combination  of  coplanar 
positions  and  precision  points. 


Motion 


CHAPTER  VII 

NINE  MULTIPLY  SEPARATED  PRECISION  POINTS 
IN  COPLANAR  MOTION 


From  Chapter  II  it  was  shown  that  five  positions  of  the 
moving  plane  determine  a maximum  of  four  Burmester  pairs. 

For  the  four  Burmester  pairsthere  are  six  combinations  of 
two  pairs.  These  six  combinations  provide  six  coupler 
curves  corresponding  to  the  five  prescribed  positions  of 
the  moving  plane.'  If  one  employs  the  five  position 
criteria  for  synthesizing  nine  precision  points  of  the 
moving  plane,  the  number  of  coupler  curves  (maximum  of 
six)  presents  somewhat  of  a problem  in  optimal  numerical 
methods.  This  problem  is  the  correspondence  of  the 
initial  set  of  mechanisms  (maximum  of  six)  to  the  next 
resolved  set  of  mechanisms  (also  a maximum  of  six) . As 
shown  in  Figure  (7-1)  for  the  first  resolved  set  1A, 

2A,  3A  and  4A  and  the  second  resolved  set  IB  and  2B, 
one  can  not  distinguish  if  link  IB  is  the  result  of  a 
transition  from  1A  or  2A  or  if  2B  is  due  to  a transition 
from  2A,  3A  or  4 A . if  the  wrong  interpretation  is  made, 
the  past  history  of  the  initial  set  of  mechanisms  becomes 
insignificant  causing  the  convergence  criteria  to  be 
unstable.  Other  problems  entail  the  measure  of  the  error; 
i.e.,  it  becomes  necessary  to  analyze  the  mechanisms' coupler 


TWO  BURME3TER  SETS 


VA 


{7  -.1)  Deviation  of  a Burmester  Set  fOr 
Different  Input 
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curves.  The  analysis  of  the  coupler  curve  function  in 
terms  of  the  error  of  the  prescribed  precision  points  is 
itself  a very  involved  problem.  Therefore,  it  is  not 
desirable  to  employ  the  five  position  technique  in 
synthesizing  nine  precision  points  of  the  moving  plane. 

In  Chapter  VI,  Theorem  I shows  that  there  is  a unique 
coupler  curve  for  six  point-positions  of  the  moving  plane. 
As  will  be  shown  in  what  follows,  this  uniqueness  allows 
the  kineoatician  to  look  at  the  nine  precision  point 
problem  in  a purely  analytical  manner  in  terms  of  five 
functions  and  five  unknows.  Although  the  following 
algorithm  will  be  given  in  a general  sense  as  made  up  of 
eight  functions  aid  eight  unknowns,  the  problem  can  be 
reduced  to  five  functions  and  five  unkUMms*  using  trigono- 
metric identities. 

The  six  point-position  problem  was  established  by 
equations  (6-10)  as  two  quadractic  functions 


(7-1) 


epq  " 

are  given  by  equations  (6-11).  Also  from  Chapter  VI,  the 
singularity  constraint  functions  for  the  quadratics  are 
given  as 


r-2) 


A.«*  “ *<“ik»  8ik>  v±) 


depending  on  the  case  study.  For  the  six  point-position 
synthesis,  the  linear  set  is  given  by  equations  (6-7) 


VZ  = 
7Z  = 
7Z  = 


VZ  = 7 Z + 7 Z 


7p,  “£(^1  Z 1 34 

For  nine  precision  points  there  must  also  be  four 
comparable  linear  sets  for  Z in  terms  of  Z 

The  first  relation  of  equations  (7-3)  defines  the  line 
passing  thru  the  Burmester  points.  Therefore,  for  the 
four  linear  sets  the  first  relation  of  each  set  must  all 


be  coincident  or  identical.  By  prescripting  7 as  the 
last  prescribed  position  of  each  of  the  four  linear 
sets,  all  having  positions  0,1, 2, 3, 4 in  common,  gives 


Vz  = 
7Z  = 
7Z  = 
7Z  = 


From  these  four  relations  there  are  six  singularity 
functions  expressed  by 
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“■■"p,  - * v 


bmi  = -and  * d 
bm,  =Snd,  + d21 
bm,  " V2(«£  ♦ 8‘)dn 
+BndI1  - d;1 


^ STS£  derivatives1 


U^iSbii 


10S 


c t = -g  d + di3 

Cmi  = ~“nd»  + d,j 
c = 1/2 (a*  + 62)d 


(7-7) 


The  coefficients  d are  expressed  in  terms  of  the  first 


five  common  positions  for  i = 1,2, 3, 4, 5 and  are  given  by 


\ \ ».J 

\ 


The  singularity 


(7-8) 

may  be  reduced  to 


FJ(An  ) 


“ -&Moil 


- cnb12>1  * (b^o 
♦ (b|jCii  - eijbii)» 


= -<J> 


= - (b  c 


by  satisfying  the  identity 
sin!y.  + cossYi  = 1. 


i =1,2,3, 4,5 

{ 7—9J 


Combining  the  singularity  functions  (7-2)  with  (7-9) 


W 5 0 

a t-  >Hji  ’ 

W - 3 


Therefore,  the  nine  precision  point  synthesis  problem  has 
been  resolved  into  a five  parameter  solution  and  allows 
the  kinematician  to  look  at  the  problem  as  a purely 
analytical  function.  Figure  (7-2)  illustrates  a mechanism 
satisfying  nine  prescribed  precision  points  for  the 


By  employing  the  trigonafetric  identity  (7-10),  the 
method  is  restricted  to  a set  of  nine  precision  points 
having  a minimum  of  three  finite  angles  or  four  finite 
positions.  The  method  may,  however,  be  extended  to  include 
all  orders  of  contacts  by  lineraly  restricting  the  absolute 
values  of  the  derivatives.  It  is  felt  that  these 


108 


103 


Previous  research  by  Primrose,  Freudenstein  and 
Sar.dor  [14  J,  Tesar  and  Sparks  [15  J and  Sparks  [16] 
investigated  the  cases  of  five  and  six  symmetrical 
displacements  of  the  moving  plane.  As  shown  by  these 
researchers,  when  the  moving  planes'  positions  have  a 
characteristic  symmetry  about  some  axis,  the  corresponding 
Burmester  pairs  are  located  in  a characteristic  pattern 
about  the  axis  of  symmetry.  Since  a complete  description 
of  symmetry  can  be  found  in  references  [15]  and  [16]  no 
presentation  of  the  reductions  of  the  conics,  equations 
(2-10),  will  be  presented  in  this  text.  This  text  will, 
however,  extend  the  requirements  for  symmetry  to  nine 
precision  points  [ten  due  to  symmetry) . 

Or.o  may  employ  the  five  position  algorithm  found  in 
Chapter  II  in  synthesising  nine  symmetrical  precision 
points  of  the  moving  plane.  As  will  be  shown  in  what 
follows,  this  is  a result  of  conditions  imposed  on  four 
of  the  positions  located  on  one  side  of  V axis.  For 
simplicity  it  is  advantageous  to  separate  the  symmetrical 
mechanisms  into  three  distinct  cases  and  their  require- 


The  Burmester  points  and  their  centers  are  located 
on  the  V axis.  These  requirements  are  given  by 


conditions 


satisfied  whet 


*1 


The  Burmester  points  and  their  centers  are  symmet- 
rically located  about  the  V axis.  These  requirements 
are  given  by 


shown  that 


conditions 


satisfied 


Ill 


The  Burmester  points  are  located  on  the  U axis 
and  their  centers  on  the  V axis.  These  requirements 
are  given  by 


It  can  be  shown  that  these  conditions  are  satisfied 


In  the  synthesis  of  nine  symmetrical  precision  points 
of  the  moving  plane  it  is  necessary  to  satisfy  the 
conditions  of  one  of  the  three  cases  for  the  angles  y.  or 
the  derivative's  magnitude.  Satisfying  these  conditions 
will  give  the  imago  displacements  in  the  opposite 
coordinates  of  the  fixed  plane  as  illustrated  in  figure 
(7-3). 
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Figure  (7-3)  Symmetrical  Displacements 
Moving  Plane 


APPENDIX 


EXPANSION  OP  THE  CIRCLEPOINT 
EQUATION  BY  DETERMINANTS 


Tile  problem  of  expanding  the  compound  square  matrix 
for  four  multiply  separated  positions  in  ooplanar  motion 
is  best  represented  by  a method  of  determinants. 

The  compound  square  matrix  for  four  positions  is 
given  by 


(A-1)J 


where  the  coefficients  A are  found  in  Table  2-1.  A resulting 
general  compound  cubic  may  bo  written  as 


(A-2) 


The  coefficients 


E are  obtained  by  expanding  the  compound 


* I* 


\J* 
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cubic  (circlepoint  equation)  becomes 


<H(U  + Hsv)  (u!  + V2)  + H u2  * H V2 


4) 


(A-5) 


asymptotes  (one 
position  cubic) 


.•quation  (A- 4 ) the  circular  property 
v2)  determines  the  number  of  real 
: real  asymptote  for  the  general  four 
for  the  function. 
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